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Abstract. We find that the modified sine-Gordon equation belonging to the class of 
^SJ , the soliton equations describes the propagation of extremely short transverse acoustic 

pulses through the low-temperature crystal containing paramagnetic impurities with 
effective spin 5 = ^ in the Voigt geometry case. The features of nonlinear dynamics 
I of strain field and effective spins, which correspond to the different kinds of acoustic 

^ ■ solitons, are studied. 
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The development of physical acoustics has led to the appearance of technical tools 
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O ' of producing and measuring acoustic pulses about 10-10^ ps in duration [U The 
\0 ■ characteristics of such pulses are very perspective for diagnostics of fast processes 
^ ' and spectroscopy of solids. This attracts large attention to theoretical study of 
O ■ the interaction of picosecond acoustic pulses with paramagnetic crystals and other 
nonlinear media [3-8]. Usually, the semiclassical approach is employed to derive 
^ ' the equations governing the evolution of acoustic pulses. Some of these equations 
■ occur to be integrable with the help of the inverse scattering transformation (1ST) 
method [9l[in]. In particular, the systems of integrable equations that generalize well- 
known integrable models of nonlinear coherent optics [11] describe the propagation of 
transverse- longitudinal picosecond pulses [6], [7] . 

The duration of picosecond acoustic pulses may be comparable with the oscillation 
period of the quantum transitions involved into the interaction. Following well-known 
parallels between the nonlinear phenomena in coherent optics and physical acoustics 
[T2| [13], one has to treat acoustic pulses in this case as extremely short pulses [HI [Ti] . 
However, it is necessary in so doing to take into account essential difference between 
acoustic and optical waves. The linear velocities of the components of the former can 
differ significantly [15]. Thus, the longitudinal component velocity is normally much 
higher than the transverse ones. The nonlinear interaction of these components is 
weak in that case, and, consequently, longitudinal and transverse picosecond acoustic 
pulses propagate independently. At the same time, transverse components can interact 
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efficiently since their linear velocities are equal under propagation along the acoustic 
symmetry axis of the crystal. 

In this paper we investigate the nonlinear dynamics of the acoustic extremely 
short pulses in the low-temperature paramagnetic crystal in the external magnetic 
field presence. In accordance with above mentioned parallels between coherent optics 
and physical acoustics, we apply here the spectral overlap approximation [16]. This 
approximation is based on condition 

e = (uoTpy < 1, (1) 

where ujq is the characteristic frequency of quantum transitions created by the external 
field; Tp is the pulse duration. The main aim of the present article is to clarify the role 
of nonlinear interaction of the acoustic pulse components. We suppose for this reason 
that the pulses are especially transverse. 

Let a tetragonal (or cubic) crystal contain paramagnetic impurities with effective 
spin S = ^. Assume that the Cartesian axes x, y and z are aligned with symmetry 
axes of the crystal. Let the transverse acoustic pulse propagate along the x axis and 
the external magnetic field B be parallel to the z axis (Voigt geometry). Consider the 
one-dimensional case with dynamical variables depending on coordinate x and time t 
only. Then, the Hamiltonian H of the spin-elastic interaction has the form [13] 

H = [cr^ + F^^Sy^cTy + F^^S^^a^] . (2) 

Here uq = gfiBB/h is the frequency of the Zeeman splitting of the Kramers doublets; g 
is the Lande factor; /xb is the Bohr magneton; B = \B\; Sy^ = duy/dx and £zx = duz/dx 
are the components of the strain tensor; Uy and Uz are the Cartesian components of the 
local displacement vector n; F44 = g~^{dgyx/d£yx)o and F55 = g^^{dgzx/d£zx)o are the 
components of the tensor of the spin-elastic interaction (in Voigt notation; subscript 
"0" means differentiation at the absence of acoustic pulse); gjk are the components of 
the Lande tensor; dy and dz are the Pauli matrices; % is the Planck constant. From the 
microscopic point of view, the spin-elastic coupling appears in the case S = \ due to 
the modulation of the Lande tensor components by the strain field [T3] . 

In order to achieve fairly efficient interaction between paramagnetic impurities and 
strain field, the Zeeman splitting energy must exceed the thermal one. This implies that 
paramagnetic crystal has to be at helium temperatures, as it was in the experiments on 
acoustic self-induced transparency p2J- In that case the self-absorption of hypersound 
with frequency 10^ GHz (or the picosecond acoustic pulses) due to anharmonicity, 
defects, etc. is appreciably lower than the acoustic absorption due to the presence 
of paramagnetic impurities [13]. Hence, the self-absorption effect playing important 
role under the room temperatures can be ignored in our case. Also, characteristic phase 
relaxation time for transitions within the Zeeman multiplets is 10~^-10~^ s, and the 
energy relaxation time is much longer under such conditions [12]. We neglect these 
dissipative effects in what follows because the duration of the pulses considered is much 
shorter than all the relaxation times. 
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According to the general scheme of the semiclassical approach, we describe the 
evolution of effective spins by the equation on density matrix p: 

i^f = [^,P]- (3) 

On the other hand, the elastic pulse field obeys the classical Hamiltonian equation for 
continuous medium: 

dp 6 



dt 5u 



n<H>dr), (4) 



^ = ^{H.+ ln<H>dr], (5) 



dt 6p 

where p is the momentum density of the local displacement of the crystal; 

H.= lJ + pa' {£l + £lA dr (6) 

is the Hamiltonian of the free strain field; p is the average density of the crystal; n is 
the concentration of paramagnetic ions; <H> = Tr(/3if) is the quantum average value 
of H] a is the linear velocity of transverse acoustic waves. The integration is carried out 
over the crystal volume. 

Let us introduce the Bloch variables 

. . ^ P21 + Pl2 y ^ P2l - Pl2 yr ^ P22 " Pll 

2 ' 2i ' 2 ' 

where pjk (j, k = 1,2) are the elements of the density matrix. Then ([3D gives 

dU 

— = [uj^ + ^;)V + nyW, (7) 

dV 

— = -{uJo + n,)U, (8) 

-n,u, (9) 



dt 



where 



— ^oF44£yx, — ^oF^b^zx- 

With (ED, (!4D-(I6D we obtain 

d'Qy ^2d'^y _ nhulFld'V 
dt"^ dx^ 4p dx^ ' 



dt"^ dx^ Ap dx"^ 

Equations flTD- ffTTp describe the interaction of the transverse strain field with the 
paramagnetic crystal in the Voigt geometry case. As it is seen from ([7D-([9D, y-component 
VLy of the acoustic pulse causes quantum transitions between the Zeeman sublevels, 
whereas z-component VL^ shifts dynamically their frequency. For transverse acoustic 
pulse propagating along the z axis (Faraday geometry), both components of the pulse 
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excite quantum transitions only. The spin-elastic interaction between the components 
leads in this case to the rotation of the polarization plane of the pulse |1] (acoustic 
Faraday effect). 
If we put 

S = W + iU, 
then (C]) and (El) yield 

— = i{ujo + Q,)V + iQyS. (12) 

Let us assume that Tp ~ 10 ps and the orders of ujq and are comparable. Taking 
uo ~ 10^0 s-^ (that is 5 ~ lO^Gs) H [121 [131, we see that condition ^ is vahd. In 
that case the first term in the rhs of ([T2|) can be neglected in the approximation of 
zeroth-order with respect to e [T7]. Then we have 

U = Wo sin e, W = Wo cos 9, (13) 



or 



where 



fnydt', (14) 

J to 

Wo {\Wo\ < 1/2) is the inversion of population of the spin sublevels in the acoustic pulse 
absence. Substitution ([T3ll into (JH]) gives 
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^ = -Wo{uJo + n,)sm9. (15) 
at 

To simplify further the equations we deal with, let us carry out some numerical 
estimations. Assuming W ~ U, d/dt ~ l/xp we find from that fly ~ l/r^. 
Therefore, the ratio rjy of the rhs of equation f[TO|) to the terms in its Ihs is estimated 
as riy ~ ^/enh,uJQF'l^/ Apa? . The value of similar parameter of ( [TTl) is estimated 
as Tjz ~ y/enhujQF^^/ Apa? . For paramagnetic ions Co^^ in cubic crystal MgO at 
helium temperatures we use the following experimental data [H [IBJ: n ~ lO^^cm"^, 
ujQ ~ 10^° s-\ p ~ 1 g/cm^, a ~ 5 ■ 10^ cm/s, and F44 ~ F55 ~ 10^. If Tp ~ IQ-^^ s, then 
Vy ^ Vz ^ 10~2. Since parameters r]y and rj^ are much less than unity, we shall reduce 
the order of derivatives in ( ITOl) and ( [TTl) with the help of the unidirectional propagation 
approximation [T8] . 

Having introduced new independent variables r = t — x/a and ( = rjx, where 
r] = max{r]y, rj^), we obtain 

d d d 1 d d 

dt 9r ' dx adr ^ d(' 
In the first order in 77, we write 

Q2 1 Q2 Q2 Q2 1 Q2 



-2 



dx"^ dr"^ a drdC, ' dx"^ a? dr'^ 
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for the Ihs and rhs of equations fllUp and fill I) , respectively. Integration of the wave 
equations obtained in this way with respect to r, substitution of expressions f|T3|) and 
taking into account ( |T5i) give us the following system in the terms of variables r and x: 
dQ 

" -Py{uJo + n,)sme, (16) 



dx 



Pzfly sin 9, 



dx ^' 

,2 T?2 /{Q^„3\ 



[17) 



where (3y = -WonUujlFl/ {9>pa^), (5, = f3yFl/ F^. 

Equations f|T6l) and f|T7j) possess the integral of motion: 



-^44 



where function /(r) is determined by the boundary conditions. The similar integral was 
revealed in [?]. Defining new variables 



1 + f{f)/ujl df, 



1 + f\r)/ul 



one can prove that /(r) is supposed equal to zero without loss of the generality [Tj. 
Then, we find from ([H 



where 



'55 



(19) 



1^0-^44 

(It is seen that inequality \VLz 
obtain 



< 2uJo is fulfilled.) Finally, using (fMl), (fT6l) and (fT9l). we 



dxdr 



\ 



l-r2 



(20) 



This equation is reduced to the famous sine-Gordon (SG) equation P, Unj if 
Tc = 0. Equation fl20l) with Tc 7^ is known as the modified SG (mSG) equation 
[19-22] and belongs to the class of equations integrable by the 1ST method. Its first 
physical application was found recently in [23], where (l20l) was shown to describe the 
propagation of electromagnetic extremely short pulses through the anisotropic media. In 
[19-22], this equation was derived in the course of mathematical study of the Backlund 
transformation of the SG equation. 

Being integrable with the help of the 1ST method, (l20ll admits the zero curvature 
representation 

dL dA 0. _ 
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r. 



A 



LJoPy ( — ircsin6' Ae^^ 

-n 

Ae sin 6 



and A is the spectral parameter. Equation (!2T]) is nothing but the compatibihty condition 
of the following Lax pair 

OX 

where ^ = ^(A, r, x) = (^i, 6)^- 

To investigate the nonlinear dynamic of the transverse strain field components and 
effective spins, we construct the soliton solutions of fl20|) . It is well known that the multi- 
soliton solutions of the integrable equation can be found using the algebraic methods. 
Here we apply the Darboux transformation (DT) technique [23]. Let (f = {ifi,(f2)'^ 
be a solution of (122|1 with A = Tp. The Lax pair (!22l) is covariant with respect to DT 
{6, 6, 0} {fi, 6, 0} of the form 

6 = (A^i - Wi6/^2) exp[i(0~ - e)/2], 
6 = (A6 - W2ei/^i) exp[i(e - ^~)/2], 

^~ = ^ + iln^4^^^. (23) 

This implies that relation (l23ll gives us new solution ^ of the mSG equation (l20l) if 6' is 
its known solution and is a solution of the Lax pair. 

In the zero background case (i.e., 6 = 0), we obtain from ( |23l) the following 
expression for the one-soliton solution of the mSG equation: 

q- tanh x 
U = Z arccos 1= , 

where q = Tc/Tp, x = {t — x / v) / Tp, A = 1 — 2gtanhx + g^. Velocity v of the soliton and 
its free parameter Tp defining the duration are connected by the relation 

v~'^ = + uJoPyTp. 

The corresponding formula for y-component of the transverse strain field is 

_ 2 sech X 1 -gtanhx . . 

For "time area" Ay = Qy dt of this component of the acoustic pulse we find 



Ay 



±2tt for \Tp\ > Tc 
for \Tp\ < Tc 
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The last formula indicates that the acoustic extremely short pulses are divided into 
two families. The family with Ay = ±27r exists for the SG equation also and corresponds 
to unipolar 27r-pulses (kinks and antikinks). The pulses of the family with Ay = are 
bipolar Ovr-pulses. Unlike the breathers of the SG equation (for them Ay = as well), 
these pulses are steady-state. The solitons of this kind were called as neutral kinks in 

m- 

In the cubic crystal, one has IF55I = IF44I. Then \Tp\ < r^. due to condition ([1]), and 
neutral kinks exist only in such a crystal. In the crystals with tetragonal symmetry, 
both types of the solitons are possible. 

Let us discuss in details the properties of the acoustic solitons. In the case 
\Tp\ > V2tc, component Vly fl2^ of the unipolar one-soliton solution has a single 
maximum, whose value is smaller than 1 / Tc. Accompanying dynamics of effective spins 
is very similar to that of the SG equation: the leading edge of the pulse of ^/-component 
inverts completely the populations of the Zeeman sublevels, and the trailing one returns 
them to the initial state. The ^-component is small as compared to ujq, and its role is 
insignificant. 

Under Tc < \rp\ < V2tc, component f2y has two symmetric peaks (see solid line in 
figure la) with the largest possible amplitude l/r^ determined by (fT8|) . The peaks are 
separated by the time interval 

2 r„ arcsinh W — 1. 

V 1 — g 

The first peak inverts the populations of the spin sublevels while the component Qz 
grows in amplitude and reaches the absolute value 2a;o in the center of the soliton 
(figures lb and Ic). The component Qz has an asymmetry on polarity: it decreases 
the transition frequency {Qz < 0) and shifts the Zeeman sublevels so that the ground 
sublevel becomes excited. On account of this, the paramagnetic ions are in the ground 
state between the peaks. (The positions of the spin sublevels 1 and 2 of the Kramers 
doublet under the pulse passage are pictured in figure Ic.) When the second peak has 
come, the 2;-component vanishes reverting the mutual position of the sublevels to the 
initial state. Finally, the second peak of y-component causes the back transitions from 
excited sublevel to the ground one. 

When \Tp\ < Tc and the interval between the peaks of neutral kink surpasses its 
duration, the dynamics of the strain fields and effective spins is similar to the second 
case described above. The only difference is that the peaks of fly are opposite in sign. 
The time interval between the peaks is 

, , , I ^ 
2 Tp arccoshWl H — . 

V - 1 

If we take duration of such a soliton to be shorter, then the peaks are brought closer 
together and the degree of excitation of the paramagnetic ions decreases (see dotted 
lines in figure 1). 

When \Tp\ — > Tc, the interval between the peaks grows indefinitely large, and y- 
component (l24l) consists of a single peak with amplitude equal to 1/tc and with absolute 




Figure 1. Profiles of the components of tlie strain field and tlie population pu for 
the one-soliton pulses with Tc < \Tp\ < \/2tc (solid lines) and \Tp\ < Tc (dotted lines) 



value of time area Ay equal to vr. This case stresses especially the role of the component 
Qz of the acoustic pulses considered. The peak of Qy inverts almost completely the 
population of the spin sublevels. This state of the effective spins is unstable in the 
absence of the strain field. But, the ^-component, whose amplitude tends to 2u!q, shifts 
the levels of the Kramers doublets in a such manner that the energy of the excited 
sublevel becomes lesser than the energy of the ground one. Owing to this, the state of 
effective spins after the passage of the ^/-component peak becomes stable. 

The form of the acoustic solitons in the case /(r) 7^ (see integral (flSl) ) tends to 
their form in the case /(r) = at t ^ 00 since the pulse velocity v differs from linear 
velocity a of the transverse waves. As it follows from the previous consideration, this 
means in particular that the amplitude of the pulses is bounded, and the component of 
the strain field parallel to the external magnetic field has the asymmetry on polarity. 
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In this paper we considered the propagation of the transverse acoustic extremely 
short pulse through paramagnetic crystal in a direction perpendicular to external 
magnetic field. It was shown that the dynamics of the strain field and effective spins 
is governed by the modified sine-Gordon equation (120|) . The soliton solutions of this 
equation reveal strong nonlinear coupling between the components of the acoustic pulse. 
As a result of this, the behaviour of paramagnetic impurities and elastic fields during 
the interaction exhibits new features. 
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